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Abstract 

We give an idea of the evolution of mathematical nonlinear geomet- 
ric optics from its foundation by Lax in 1957, and present applications 
in various fields of mathematics and physics. 

1 Introduction 

Geometric optics goes back at least to the XVIIth Century, with Fermat, 
Snell and Descartes, who described the "paths" (rays) followed by the light. 
Nowadays, Physics tells us that we may reasonably replace the waves from 
Quantum Mechanics with classical particles, in the semi-classical approx- 
imation (when considering Planck's constant h, or the wavelength, as in- 
finitely small). The mathematical transcription of these problems consists in 
studying the asymptotic behavior of solutions to partial differential equations 
where different scales (represented by small parameters) are present, often in 
a high frequency oscillatory context. 

We present the first historical results of the field, and then review some 
extensions and applications of the method. We shall see how geometric 
optics applies to Maxwell's equations (from optics, ferromagnetism, . . . ), to 
the wave or Klein-Gordon equation, to fluid dynamics and plasma physics, 
to general hyperbolic systems and conservation laws, as well as to nonlinear 
Schrodinger equations, among others. Furthermore, we shall see that it does 
not apply to oscillatory problems only, but also to boundary layers, shocks 
and long waves problems. 
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Quoting some reyiews and introductory texts may be useful. A reyiew 
on nonlinear geometric optics in 1998 is due to Joly, Metiyier and Ranch 
|.TMR,99j . and the onhne book of Ranch |R,au96j is a nice introduction to 
the subject. Majda |Maj84| made major contributions in the 80ies. Good 
Physics textbooks on nonlinear optics are Newell and Moloney [NM92j and 
Boyd |Boy03| . Kalyakin |Kal89j and Hunter |Hun95j reyiew many questions, 
methods and applications about these multi-scale problems, and Whitham 
Whi99j had pioneering contributions in the modeling of nonlinear wayes. 
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2 First steps 



2.1 Linear geometric optics 

The first rigorous result in mathematical geometric optics is due to Lax 
Lax57j . who shows that strictly hyperbolic systems admit WKB (for Wentzel, 
Kramers and Brillouin) solutions. These have the phase- amplitude represen- 
tation 

The wavelength corresponds to the small parameter e > 0, and the amplitude 
admits an asymptotic (in general not convergent) expansion, 

a*^ ~ ^£:"'a„(t, x), as e — > 0. 

n>0 

In particular, Lax uses such solutions to study the well-posedness of the 
Cauchy problem. 

The geometric information is contained in the eikonal equation (of Hamilton- 
Jacobi type) determining the phase ip (plane waves correspond to phases 
which are linear in (t, x); spherical waves, to phases which are functions of t 
and only). The amplitudes a„ are solutions to transport equations along 
the rays associated with the eikonal equation. The asymptotic expansion is 
valid for times before rays focus. 

2.2 Weakly nonlinear geometric optics 

Trying to generalize Lax's approach to nonlinear systems, one immediately 
faces several problems. First, a family (a(a;)e*'^^^-*/^)^^Q is not bounded in 
any Sobolev space if* with s > 0, so that the time of existence of the 
considered exact solution may shrink to zero as e goes to zero. In the 
same spirit, the size of must be adapted, large enough to allow nonlinear 
features at first order, but sufficiently small to prevent blow-up. Finally, 
the main interest and difficulty of nonlinear models comes from interactions: 
one hopes to incorporate in the asymptotic description the self-interaction 
of a wave (leading for example to generation of harmonics) as well as the 
possibility of (resonant) interaction of several waves. 

Thus, the extension to nonlinear systems goes through profiles Un{t, x, 6) 
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which are periodic or almost periodic w.r.t. 6 G M"^, 

(1) ~ x) + e"" ^ e"f/„ {t,x,0{t,x)/e) , 

n>0 

where m is a given groundstate, and = {ipi, . . . , (pq) is a collection of phases. 

The usual strategy of nonlinear geometric optics consists in: 1-defining a 
formal solution, i.e. solvable equations for the profiles; truncating the series 
in (P) defines an approximate solution (a function approximately solution to 
the equation); 2-showing that for any initial data close to the initial value of 
the approximate solution, an exact solution exists on a time interval indepen- 
dent of e; 3{stabiUty)-showmg that the exact solution is well approximated 
by the formal one. 

In order to observe some nonlinear behavior, the magnitude of oscilla- 
tions is chosen so that cumulated effects of nonlinearities become of the same 
order as the wave on the typical time T of propagation. A fixed T (w.r.t. 
e) leads to nonlinear geometric optics. Once rescaled, the system of partial 
differential equations takes the form 

(2) L{t, X, u^ ed)u' = F{t, x, u"), 

where the operator L is in general a first order symmetric hyperbolic system 
on some domain in Mj^"^, 

d 

L(t, X, u, ed) u = edtu + Aj{t, x, u)edxM'^ + i^o(^) u 

(3) j=i 

= : Li(t, X, u, ed) u + Lo(t, x) u. 

The smooth functions F, B and Aj take values in C^, in the space of x 
N matrices, and in the space of N x N symmetric matrices, respectively. 
Furthermore, we will distinguish the hyperbolic case, when Lq = 0, and 
the dispersive case, when Lq is skew-symmetric (the dissipative case, when 
Lq + Lq > 0, will be considered only in Section [7j). This coefficient Lq reflects 
the interaction between the wave and the material medium, so that, in the 
dispersive case, the group velocity (see Q) depends on the frequency. The 
abovementioned exponent m is then defined as follows. Let J > 2 be the 
order of nonlinearities, 

\a\<J-2 ^ d:^A,,=0, \/3\<J-l ^i^u.o = 0. 
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The standard amplitude of weakly nonlinear geometric optics corresponds to 



(4) 



m 



J 



(so that m = 1 for quadratic semi- or quasilinear systems). 

As an example, Euler equations for compressible and entropic gas dy- 
namics form a quasilinear hyperbolic system (Lq = 0), and Maxwell-Bloch 
equations, a semilinear dispersive one (describing the propagation of an elec- 
tromagnetic wave [E, B) in a medium with polarization vector P; for a two- 
level quantum medium, N is the difference of the populations of the energy 
levels), which reads 



(5) 



< 



(dtB- 
dtE- 
dtN 



curl E = 
curls = 
- E ■ dtP, 



0, 



-dtP, 



^ e^d^P + P={No- N)E. 
It may be written in the form (0) for the unknowns B, E, N, P, edtP. 



2.3 Profile equations 

The formal derivation of profile equations is similar in any dimension, and 
for any number of phases. Plug expansion into equation Q, and let the 
expansion of Lu — F{u) vanish. This gives an infinite set of equations, 

L{d^-de) Uo = 0, 

and for n > 0, 

(6) L{d0-de) Un+i + Li{d) t/„ + F„(f/o,...,[/„) = 0, 

with d0- de = Ylk=i ^fbde^ ■ 

The analysis is based on formal Fourier series in 6, so that for U = 

= L{d0 ■de)U = Y^ ^(^d(a ■ 0)) f/^e^"'^ = UU = J2 TTaf/^e^^'^ 
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where tTq is the orthogonal projector on the kernel of L{id{a ■ (f)). In order to 
get non-trivial solutions, this projector is assumed not to vanish at least for 
one multi-index a. This means that the function a ■ (p satisfies the equation 

(7) detL(id(a ■ (^)) = 0, 

so that a ■ (p must satisfy one of the eikonal equations, 

dt{a ■ 0) + \k(t, X, d^{a ■ (p)) = 0. 

Here, we denote by Xk{t, x,^),l < k < N, the N eigenvalues of the symmetric 
matrix 

1 1 
-L{t, X, 0, (0, tO) = J2 ^' 0) + 

' i=i ' 

In particular, a single initial phase ip^{x) may generate eikonal phases 
(pk{t,x) (except when initial profiles are prepared, i.e. polarized on a single 
mode). 

When the linear combination a -(pis not trivial, i. e. contains linearly inde- 
pendent v^fc's, equation ((7j) expresses a resonance between these phases. Such 
resonances must involve at least three phases, and thus occur for systems of 
size at least 3. 

Profile equations are then obtained recursively, splitting (jH)) thanks to the 
projectors 11 and (1 — 11), 

(8a) (1 - n) f/o = 0, and for n > 0, 

(8b) UL^{d)UUn + IlFn = 0, 

(8c) (1-n) Un+i = -L{dip-de)-'[Li{d) Un + En]. 

Equation ()8a|) is a polarization constraint on the first profile Uq, and ()8b|l 
gives the evolution for the polarized part 11 f/„ of Un in the domain Q C M^+'^ 
considered. The first equation ()8b|) (for the first profile Uq) is nonlinear, while 
the others (for n > 1) are linear. The operator L(d<^-9e)~^ is formally defined 
by the symbol L(id(a ■ <^))~^, the partial inverse of L(zd(a ■ 0)) on its range. 

The description of the operator 11 Li 11 depends on the geometry of the 
characteristic variety, 

C = {{t,x,T,0 eQ X R^+'^ I det L{t, x, r, = 0}, 
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viewed as a different iable manifold. In the hyperbolic case (Lq = 0), this 
characteristic variety is homogeneous in (t, ^): for the wave equation, it is 
the light cone C = fix {r^ = In the dispersive case (Lq = Lg 7^ 0), only 

few harmonics of an eikonal (r, ^) are eikonal in general: for Klein-Gordon 
equation, C = x {r^ = + 1}. 

When for all {t,x) G Q, the vector d{a ■ 0){t,x) is a smooth point of a 
sheet T = X{t,^) of C, the principal part of naLi{d)TTa is simply the (scalar) 
transport operator at the group velocity, 

dX 

(9) Vait,x) = - — {t,x,d^{a ■ (fi){t,x)); 

(for the wave equation, (9g(±|^|) = ±^/|^| does not depend on |^|, whereas for 
Klein-Gordon equation, d^{±^\E,\'^ + 1) = / Vl'CP + 1 does). Precisely 
r |Lax57j . |DR97aj ). 



■KaLi{d)-na = dt + Va-d^ + 7r„(Li((9)7r„). 

This reduction clearly explains the "light rays" picture; it is however 
not necessary in the justification of the asymptotics as an example, 
Joly, Metivier and Ranch perform this justification JMR94 J in the case of 
an eigenvalue A which changes multiplicity (crossing in the characteristic 
variety), which gives a model for conical refraction of light. Moreover, this 
description is clarified by Lannes |Lan98j . in the case of isolated singular 
points of a dispersive C: the characteristic variety of TraLi{d)7ra is then the 
tangent cone to C. A systematic approach to such "algebraic lemmas" is 
given by Texier |lbx04j via perturbation theory of complex eigenvalues under 
a slightly different smoothness assumption (found by Metivier and Ranch 
MH.n.Sj to be actually the same in the hyperbolic case). This approach also 



generalizes the one of Lannes, and applies to "all order algebraic lemmas" 
(see Section EI). 



2.4 Rigorous results 

One easily suspects, in view of the term L{d0- de)~^ in (jHcj) . that resonances 
may cause troubles. That is the reason why monophase expansions (with 
profiles polarized on a single mode) have been studied first. 

In |GB69j . Choquet-Bruhat constructs formal such expansions for quasi- 
linear systems, and applies these to computations on the model of relativistic 
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perfect fluids. Single-phase weakly nonlinear geometric optics is justified in 
space dimension d, for semilinear systems, by Joly and Ranch |J R92j in the 
non-dispersive case, by Donnat and Ranch |DR97aj . |DR97bj in the disper- 
sive case, and by Gues |Gue93j for quasilinear non- dispersive systems. Their 
strategy consists in constructing an approximate solution w^pp, sum of a 
high number of smooth profiles (or with infinitely accurate asymptotics, via 
Borel's summation technique), so as to absorb the oscillations of the resid- 
ual Lulpp — F(m|pp). Existence of exact solutions and stability are then 
obtained via energy estimates on — w^pp. 

The formal study of propagation and interactions of oscillations for quasi- 
linear systems, particularly for fluid dynamics, is done in the works of Hunter 
and Keller |HK83j in the non interacting case, of Majda and Rosales [MR84j 
in the one dimensional resonant case, of Hunter, Majda and Rosales [HM R86] 
and Pego Peg88| in the multidimensional resonant case. In these works, a 



finiteness assumption is made on the number of (directions of) eikonal phases. 
But in [JR91 , Joly and Ranch show that resonances may generate dense os- 
cillations in the characteristic variety. 

This finiteness hypothesis is shown unnecessary in |,IMR98bj . where Joly, 
Metivier and Ranch justify rigorously the weakly nonlinear geometric optics 
approximation in space dimension 1, for both semi- and quasilinear systems 
(before shock formation), under a weaker transversality assumption between 
real combinations s of phases and the propagation fields X^. {X^s = 0, or 
XfcS 7^ almost everywhere on the domain Q). They use the method initi- 
ated by Joly |Jol83j for semilinear multidimensional systems with constant 
coefficients and linear phases, estimating (by stationary and non-stationary 
phase method) the difference — e"^UQ{t,x,0/e) between the fc-th Picard 
iterates obtained from the resolution of Q and of the first profile equation, 
respectively. This leads to a first order approximation, 

u%t, x) = e'^Uoit, X, 0{t, x)/e) + 0(5"^), 

which is valid in V for all p < 00. It holds in L°° when transversality 
is imposed (almost everywhere) along integral curves of each propagation 
field Xfc. Cues shows |Gue95aj that on the contrary, when weak resonances 
are present (a linear combination of phases is eikonal on a set with one 
dimensional positive measure, but is a.e. not eikonal on fi), they may lead 
to the creation of stationary waves preventing from L°° approximation. 

Concerning results in space dimension one, we have the earlier works (for 
linear phases) of Kalyakin jKal89j , and also of Tartar et al. [Iar81j , |l'ar84j , 
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|MPT85j . by the compensated compactness method apphed to Young mea- 
sures -hnks between nonhnear geometric optics and compensated compact- 
ness are very clearly enlighted in |MS98j and |.IMR95cj . 

Furthermore, still in space dimension one, several papers deal with (global) 
weak solutions (with bounded variations) of conservation laws. The formal 
study is due to Majda and Rosales jMR84j : DiPerna and Majda obtain LfL]. 
global in time asymptotics -even when shocks occur-, using L} stability and 
BV decay properties of solutions (constructed by Glimm's scheme) to con- 
servation laws, for BV initial data with compact support. This result is 
extended by Schochet Sch94bJ to periodic BV data using WKB expansions, 
and then by Cheverry |Che96aj , |Che97j to general B V data. Sable- Tougeron 
|ST96j treats the initial-boundary value problem, for linear non-interacting 
phases, and data with compact support. 

Now, the multidimensional case for semi- or quasilinear interacting waves 
is treated by Joly, Metivier and Ranch [JMR95a]. In space dimension greater 
or equal to 2, the new phenomenon is the focusing of oscillations, correspond- 
ing to singularities of phases (see Section E)). It may lead to blow-up, and to 
the ill-posedness of the Cauchy problem. Joly, Metivier and Ranch distin- 
guish "direct focusing" of principal phases and "hidden focusing" resulting 
from several interactions. They define a coherence criterium so as to separate 
clearly eikonal and non-eikonal phases. 

Definition 2.1. A real vector space $ C C°°{VL) is L-coherent when for all 
ip G $\{0}, one of the following conditions holds: 

(i) det L(t, X, dipit^ x)) = and dip{t^ x) 7^ for every {t, x) G Q, 

(ii) det L{t,x, dp (t,x)) does not vanish on Q. 

The typical example of L-coherence is for linear phases and constant 
coefficient operator L. Coherence is not easy to check on a given set of 
phases, but it allows to prove (locally, on truncated cones Qt = {(^,2;) | < 
t < T , \x\+t/6 < p}, for small S, p > 0), for a symmetric hyperbolic system 
with constant multiplicity, the following result. 

Theorem 2.1. Consider initial phases p^ = {p^, . . . ,p^) such that dpi 
does not vanish on Qq. Consider eikonal phases p = {pi, . . . , pq) generated 
by these initial phases (V/c, pk{0,x) = p1{x)), and assume that their linear 
span $ is L-coherent, and contains a timelike phase Pq: po{0,x) = 0, and 
dtPo does not vanish. Then, for any initial data 

(10) h%x) =e'''n%x,pye), 



9 



with (?i'^)o<e<i (I bounded family in the Sobolev space H^{Qq x T'') for some 
s > 1 + {d + q)/2, we have 

(i) There is T > such that for all e g]0, 1], the Cauchy problem (0), (fTUIl 
has a unique solution G C^{Qt)- 

(ii) There is a profile U{t, x, r, 6) G C^{Qt x M x T'^), almost periodic in (r, 6) 
(and determined by the profile equations (jSaj), (jsHj) ) such that 

u%t, x) - e'^Uit, X, <fo/e, 0/e) = o(e'") in 

The new idea in their proof consists in looking for exact solutions of the 
form 

u'{t,x) = eV{t,x,if/e). 

It is then sufficient (for to be a solution to 0) that the profile be 
solution to the singular system 

L,{dt,^)W + -^ i^LM^k)de,+L^ W = F'{t,xM'), 

and L-coherence precisely allows energy estimates in H'^ . 

The result of Joly, Metivier and Ranch on singular systems generalizes 
the ones of Klainerman and Majda |KM81j . |KM82j and Schochet [Sch94a] . 



It is worth noting that the analysis of such singular systems appears in many 
multi-scale problems. In particular, for rapidly rotating fiuids in oceanol- 
ogy, meteorology, and Magneto- Hydro-Dynamics; see [BMN96j, |Gre97bj . 
Gal98j ■ |MS01j : several applications are described by Klein [ KleOOj . This 



method also provides an efficient tool for numerics: see the works of Colin 
and Nkonga about propagation of wavetrains |CN04bj and pulses (see sec- 



tion ini) |CJN05j in optical media, and Colin and Torri |CT05j about pulse 
propagation over diffractive scales (see section ing . 

A work close to the method exposed here may be found in the paper 
|CCarj by Colin and Colin, studying rigorously and numerically Raman scat- 
tering for a semi-classical quasilinear Zhakarov system from plasma physics 
in the weakly nonlinear geometric optics regime, with three-wave resonances. 
There, the system lacks hyperbolicity, which is compensated by dispersion. 

It is in fact possible to achieve infinitely accurate asymptotics for these 
multidimensional nonlinear interacting waves with nonlinear phases, by adding 
to the coherence assumption on phases a generically satisfied "not too small 
divisors" hypothesis (or Diophantine hypothesis on wave vectors, in the case 
of linear phases) , see |JMR93a] . 
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3 Other kinds of profiles 



As we have seen, the natural dependence of the profiles U(t, x, 0) on the fast 
variable ^ G M*^ is periodic, quasi- or almost periodic. This allows to define an 
oscillating spectrum |,Tol88j . |.TMR,94j which is localized on the frequencies 
given by the Fourier transform of U{t,x,-). Now, WKB aymptotics with 
profiles having other behaviors w.r.t. 6 may be relevant, depending on the 
context. The formal computations are in general very similar to the ones 
of usual geometric optics, but at least, the functional tools (such as average 
operators) must be re-defined, and interactions may take a different form. 

In order to model the optical Raman scattering, for which light is emit- 
ted in a continuum of directions, Lannes |Lan98j introduces profiles with 
continuous spectra, and gives the analogue of usual rigorous weakly nonlin- 
ear geometric optics in this context, with a precise analysis of resonances. 

In fact, this formalism also includes (see Barrailh and Lannes |BL02j ) the 
one of "ultrashort pulses" (from laser physics) considered by Alterman and 
Ranch |AR02j . where profiles have a compact support in 9 (see also sectionsEl 
and ini about focusing and diffraction of pulses). 

Sometimes, profiles may also have different limits as the variable 6 E M. 
goes to ±oo, in order to match boundary conditions (see Section 13^ . or to 
describe transitions like the ones of |Dum04j . between light and shadow (see 
Sectional); this is also the case of the solitary waves in |Gue95aj . 



4 Caustics 

In space dimension greater or equal to 2, singularities appear, even in the 
case of linear geometric optics, in the resolution of the eikonal equation: when 
rays have an envelope (called a caustic; the example of the cusp {(t,xi) e 
[0, oo[x]R I t^/^ = x1^^ + 1} from the wave equation in space dimension 2 and 
(p{0,x) = X2 + xf is shown on Figure HJ, the second derivative of the phase 
{p becomes singular. In this case, the amplitude is also singular |Lud66j . 



In order to solve the eikonal equation globally in time, one considers 
the Lagrangian manifold folliated by the Hamiltonian flow associated with 
L(t, x,T,^), starting from points (0, x, 0, dip{0, x)). Caustic crossing induces 
a phase shift determined by Maslov's index, and generates new phases (see 
Duistermaat |Dui74j : one phase before the cusp of Figure ^ corresponds to 
three phases beyond the cusp). 
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Following these ideas, Joly, Metivier and Ranch have stndied the canstic 
crossing for semilinear geometric optics. In |JMR95b] . they show that, for 
superlinear nonlinearities and focusing at a point, oscillations may lead to 
explosion, whereas in the case of dissipative equations (for which exact solu- 
tions are globally defined), they may be absorbed {i.e. only a non-oscillating 
term remains after reflection on the caustic). In JMROOa , for dissipative 
equations, they extend this result to general caustics. They exhibit a critical 
exponent pc, defined by the geometry of the caustic so that, if the nonlinear- 
ity is stronger than lul^" at infinity, absorption occurs, and else, oscillations 
persist. In [JMR96 , for uniformly Lipschitzean nonlinearities (for which, 
again, exact solutions are globally defined), they show that oscillations go 
through the caustic. 




Figure 1: The cusp. Through each {t,x) point pass only one ray before the 
caustic, and three rays beyond the caustic. 

More precisely, to a single-phase expansion before the caustic (case con- 
sidered here just for notational simplicity; several initial phases are allowed), 

u'{t, x) = e™W(t, X, if/e) + oil), 
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corresponds a multiphase expansion beyond the caustic, 

j 

with phases (pj defined as in the hnear case, matching ip on the caustic. In 
particular, this provides examples of asymptotics with non-coherent phases. 
This approximation is in in |JMR95c] and |JMR96j . and in all L^'s, p < oo, 
in [JMROOaj . It is based on the oscillating integral representation 

u'it, x) = e^'uit, x) + y e*^(*'^'?^'«)/^^(t, x, ip/e) dy d^ + o(l) 

= : e"'u{t,x)+e'''r{A) + o{l). 

Outside the caustic, expanding a smooth A(t, x, 6) in Fourier series in 6 and 
performing the usual stationary phase asymptotics (the assumption that the 
phases tpj are non-resonant is needed) gives 

3 

where the sum corresponds to the several pre-images of critical points of 
the phase, and Jj{A){t,x) is the evaluation at {t,x,ipj/e) of some profile 
Uj{t,x,6) obtained from A via a Hilbert transform. 

Here, the amplitude A belongs to L^, so that the profiles Uj are in fact 
better defined as "weak profiles" (in particular, Uj{t,x,(pj/e) has an asymp- 
totic sense only). They are weak limits of m^, using Jj{B) as test function, 
for smooth B(t, x, 9), 

\fB{t,x,e), J WJ]{B)dx^ J UjBdxdd. 

Profile equations are then obtained as weak limits of the original equation, 
using a nice "no propagated oscillations" lemma f |JMR00aj . section 5.2) for 
nonlinearities. Strong convergence then follows from energy estimates and 
U estimates on oscillatory integrals. 

These results are refined (in L°° instead of L^) and extended to conserva- 
tive systems by Carles, for simpler geometries. Essentially, Carles considers 
focusing on a point for the wave equation, 

(11) Uu' + a\dtu'\P-^dtu' = p>l, aGC, 
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{a > 0, a < 0, a E iM. corresponds to the dissipative, accretive and conserva- 
tive case, respectively) with radially symmetric initial data, or the semiclas- 
sical nonlinear Schrodinger equation (NLS), 

(12) iedtu' + ^e^Au' = e''\u'\'^''u', a > 1, a > 0, 

with initial data of the form 

Hunter and Keller |HK87j give a formal classification of the qualitative 
properties of weakly nonlinear geometric optics for (jllj) . separating linear 
and nonlinear propagation, and linear and nonlinear effects of the caustic. 
Carles has rigorously explained this classification. In |Car98j . for 1 < p < 2, 
he improves the results of j-IMRnOaj r |.IMR95bj corresponds to a > 0, p > 2) 
with an L°° description of radial waves in M^. He shows that the profiles are 
really singular, and that new amplitudes (of size e^~^) must be added to the 
one in |HK87j near the caustic. 

In a series of papers |(yRn2j . |(yR04aj . |CR,04bj . Carles and Ranch study 
the focusing of spherical pulse solutions to (fTT|) in space dimension 3, for 
various powers p and sizes (e™) of initial data. In particular, they show how 
pulses get wider after refiection on the caustic; they also treat the analogue 
of |JMR95b] (a e M) and find the same results (dissipation or growth) as for 
wavetrains. Carles and Lannes generalize these results |CL03aj to dispersive 
semilinear systems such as Klein-Gordon and Hartree's equations. Again for 
ffTTj) . they obtain [CLOSbj the focusing (at t = 1) of a pulse with "random 
phase shift" ln((l — t)/e) -in the conservative case, for nonlinear propagation 
and a nonlinear caustic. 

Concerning NLS (jl2j) . Carles investigates all the behaviors with different 
powers a > 1 and a > 0. In jCarOlj . he shows the possibility of "random 



phase shift". In [CarOObj {a = da > 1, in space dimension d), he notices 
that nonlinear effects take place near the focal point only; thus, in |CFKG03] , 
|CKarj . he answers with his co-authors the converse question "for which kind 
of initial data can we get nonlinear effects at first order?" . He also considers 
|CarOOaj focusing on a line instead of a point. 

Other extensions of [CarOObj concern the description of Bose-Einstein 



condensates, via the addition to (fT^ of a polynomial potential |Car03j . 
|C;M()4j . |C;N()4aj . or Hartree-type nonlinearities ((l/|x|) \u\'^)u |C;MSarj 
for the Schrodinger-Poisson model. 
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The method of Carles consists in a precise control of the solution outside 
the caustic via WKB expansions and Gagliardo-Nirenberg estimates, then 
rescaling the problem around the focal point so as to use nonlinear scattering 
arguments (with short-range or long-range interpretation, according to the 
power of the nonlinearity) . 

We may mention that caustics are subject to many numerical works, such 
as the ones of Benamou et al, [FSOO], |BH()2j . |Ben():-ij . |BLSS():-ij . |BLSS()4j . 
based on a Eulerian -or ray tracing- approach. See also the review of Gosse 
|Gos05j -about (NLS), mostly with kinetic methods- and references therein. 



5 Boundary problems 
5.1 Reflection on a boundary 

The formal study of reflections of weakly nonlinear waves is done by Artola 
and Majda in jMA88j . In the case of transverse reflection, rigorous results 
are due to Chikhi |Chi91j for equations with two speeds of propagation, and 
to Williams |Wil93] for general resonant interactions. 

The first order tangential reflection corresponds to diffractive points in 
the cotangent bundle of the boundary, and leads to the formation of a shadow 
zone (in fact, a caustic forms in the obstacle, so that phases become singu- 
lar at the tangency point). The linear case is treated by Friedlander |Fri76j 
and Melrose |Mel75aj . |Mel 75bJ. Weakly nonlinear geometric optics is justi- 



fied, at first order, with nonlinear incident and reflected phases, by Cheverry 
|Che96bj for globally Lipschitz nonlinearities (constructing smooth profiles, 
solutions to singular ODE's along the broken flow of the rays), and by Du- 
mas |Dum02j for a dissipative nonlinear Klein-Gordon equation (using weak 
profiles as described in section HJ. 

Higher order tangency is studied by Williams |Wil93j . |WilOOj . in the 



case of a simpler geometry (flat boundary, constant coefficient operators and 
linear phases). He obtains L°° asymptotics at all orders, incorporating in the 
description boundary layer profiles (see section 15. 2|) involving a third scale 
y/e: profile equations are then of Schrodinger type (see section IH)). 
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5.2 Boundary layers 

Boundary layers are a usual feature in the study on partial differential equa- 
tions on domains with boundaries: they stem from large variations of the 
solution in a small region near the boundary. They often appear in vanish- 
ing viscosity models, since the boundary conditions are not the same for the 
viscous (parabolic) and inviscid (hyperbolic) equations. 

Gues [Gue95b^ describes such boundary layers for viscous perturbations 
{—e'^S, with S{t,x,d) a second order elliptic operator) of semilinear mixed 
problems (associated with a linear hyperbolic operator L{t,x,d)). He con- 
structs exact solutions with infinitely accurate WKB asymptotics of the form 

u^{t,x) ~ ^i/e " W„ {t,x,ip/^/e,ip/e) , 

n>0 

Un (t, X, 9, Z) = U^(t, X) + Vn(t, X, 9) + Wn(t, X, z) + Xn{t, X, 9, z) , 

where the phase (p is transverse to the boundary, and functions of ^, ^ > de- 
cay rapidly at +00 (so that = lim Un)- Profile equations are analogue 

6,z— >+oo 

to the ones of geometric optics (but of parabolic type), with polarization 
conditions. When the boundary is non-characteristic for the hyperbolic op- 
erator L, the dependence on 6* = (p/y/e is not necesseray, and the boundary 
layer has size e; in the characteristic case, the boundary layer is wider, of 
size ^/e. Grenier j CilreQTaj considers the quasilinear case with characteristic 
boundary, and Grenier and Gues [GG98) . the non-characteristic quasilinear 
case. Extensions may be found in the work of Sueur for semilinear systems 
jSuearaj, |Suearcj . |Sue05bj . or quasilinear systems [Suearbj . |Sue05aj . 



When oscillations are present in the interior domain, they may interact 
with boundary layers. This occurs for example in the study of glancing oscil- 
lations by Williams jWil OOj . It also occurs for rotating fluid models (see the 
books of Greenspan IHieSOj and Pedlovsky |Ped79j about oceanography and 
meteorology, and Dormy |Dor97j about magneto- hydro-dynamics), which in 
the incompressible case typically have the form 

dtu' + ■ Vu" + ^Ru^ - £Au" = -V/, 
e 

where ]f is the pressure, is the divergence free velocity field, and the 
operator R is skew-symmetric (for example, Ru = e x u with a fixed vector 
e for Coriolis effect). Initial data are well-prepared when they correspond 
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to admissible data for the limit problem. Otherwise, fast time oscillations 
appear in the interior domain, and are taken into account in the profiles 
through a dependence on the variable t/e. For such studies in the well- 
prepared case, see jGM97j : in the ill-prepared case, see [MasOOj . |CDGG02] . 
|BDGV04] . jGV03bj; Gerard- Varet emphasizes the link with geometric optics 

in [HVna^ and [HVosI. 

This technique is also used for micromagnetism. Carbou, Fabrie and 
Gues |GFG02j describe the magnetization of a ferromagnetic medium, given 
by Landau- Lifshitz equation |LL69j . as the exchange coefficient goes to zero; 
Sanchez |San02j studies the diffraction of an electromagnetic wave by a con- 
ducting domain, as the conductivity goes to infinity. 



5.3 Shocks 

Shock waves entering the framework of weakly nonlinear geometric optics are 
weak shocks, since their amplitude goes to zero with e. Such an oscillatory 
shock is a perturbation of a non-oscillating shock and is given (in 

space dimension with x = {x',Xd)) by a hypersurface. 



n>0 



and on each side Xd ^ x'), a solution u%. to a system of conservation 
laws. 



n>0 



together with the usual Rankine-Hugoniot jump conditions. Of course, the 
unknowns are and ip'^ , so that this is a free-boundary problem. 

A formal description is due to Artola and Majda |MA88j . Rigorous results 
in space dimension one are given by Corli |Gor95j in a multiphase context, 
with curved phases, using the same "Picard iterates" method as in |JMR93b] 



(and the same transversality assumptions on phases). Williams |Wil99j ob- 
tains, in the spirit of |JMR93aj . multidimensional infinitely accurate asymp- 
totics with smooth profiles, around a planar shock {v^ = est, ifj = at), which 
is uniformly stable in the sense of Kreiss. This is done for linear phases (with 
a Diophantine assumption on their gradients), whose restrictions to the shock 
surface are entire multiples of a single phase cp^. 
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We also mention the work of Gues, Metivier, Williams and Zumbrun on 
the existence and stability of multidimensional viscous shocks |GMWZ04aj . 
These (non-oscillating) shocks are not weak ones (their size does not depend 
on a small parameter), and the system of conservation laws is perturbed by a 
viscous term — eAw^, so that the solution is smooth. The question is "are 
the shock solutions of the inviscid system approximated by solutions of the 
viscous system?". The answer is "yes" |GMWZ04bj . under the assumption 
of the existence of a spectrally stable wave profile U{x', {xd — ilj)/e) solution 
to the viscous system such that ?7(x', z) — > u^. These authors also study 

2— >lt0O 

GMWZOSbj . |GMWZ05aj the long time stability of planar shocks, using 
WKB asymptotics for the viscous problem. 

6 Three (and more) scale expansions: difFrac- 
tive optics 

So far, we have essentially presented problems where only two scales appear 
(at least at first sight). But more scales may be present (remember how a 
third scale appears with boundary layers in section l^72|) . In this section, we 
shall see that a third scale leads to supplement the transport equations (with 
finite propagation speed) for the profiles with a Schrodinger-type equation 
(with infinite propagation speed!) taking diffraction into account; this is the 
usual "paraxial approximation" of laser physics in Kerr media. 

6.1 Some three-scale problems 

In Donnat's thesis |Don94j . a model of "light cigars" is described for laser 
propagation. These are modulated oscillating wavetrains with frequency 
with anisotropic "support", i.e. of typical size 1 along a direction, and y/e 
along another. Such profiles are used by Boucheres, Colin, Nkonga, Tex- 
ier and Bourgeade |BCN"'"04] to describe the interaction of a laser with a 
gas, giving an approximation of Maxwell-Bloch equations (0) by a nonlinear 
Schrodinger-Bloch system. 

"Singular rays" of weakly nonlinear geometric optics may also be analysed 
with three-scale asymptotics. They are hypersurfaces (defined as '?/'(t, x) = 
for an unknown function tp) constituted by rays meeting an obstacle tangen- 
tially or at a corner. They model transitions between light and shadow (or 
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sound and silence, in the case of acoustic waves). Hunter |Hun88j introduces 
formal WKB expansions, 

(13) M^(t,x) ~£"^v^"f/„(t,a;,V^/v^,<^/£), 

n>0 

with profiles periodic w.r.t. 9 = ^p/e, and with limits a.s 1] = ip/^^ ±oo 
(such a may be understood as an smooth approximation of a contact dis- 
continuity, i.e. a characteristic shock along ip = 0). With an approach close 
to the one of |D.TMR,96] . Dumas |Dum04j validates this asymptotics rigor- 
ously, introducing new coherence assumptions on the phases {il),if), where 
ip = {t, ip) -it turns out that an intermediate phase ipQ = t is in general nec- 
essary for the profile equations to be well-posed. Profile equations are then 
transport equations in variables (t, x), and Schrodinger equations with time 
r = t/v^ and space variables x. 

Three-scale asymptotics of the form (|13p also describe oscillations ac- 
cording to perturbed phases (in a more restricted sense than in |l(Tre98j ). For 
initial data such as 

ul^^ = e^U%x, ^l/e) with = + V~e^\ 

the solution to the Cauchy problem cannot in general be approximated in 
L°° by a two-scale geometric optics description, but profiles as above, satis- 
fying a NLS equation on a torus, achieve such an approximation -see Dumas 

piimnabj. 

Now, we turn to the problem of long-time validity for geometric optics, 
in which the third scale is not present initially, but naturally comes up. This 
is the context where first rigorous three-scale asymptotics were proven in 
nonlinear optics. 

6.2 Long-time behavior 

A priori., asymptotics from section|21are valid only on some fixed time interval 
[0,to]- In fact, when the profiles are globally defined, uniformly bounded in 
space-time, and have H'^ norms with polynomial growth in time, Lannes 
and Ranch prove |LR01j that weakly nonlinear geometric optics is valid up 
to times t ~ Chi{l/e) -in the linear case, see the study of propagation on 
this so-called Ehrenfest time for semi-classical Schrodinger equations by De 
Bievre and Robert |DBH,():-ij . 
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Diffractive optics. However, in order to describe oscillatory waves with 
frequency 1/e on larger propagation scales, one must adapt the geometric 
optics approach (even in the linear case). First, the cumulated effects of 
nonlinearities over longer times leads to consider smaller amplitudes. The 
first natural long-time scale is then the diffractive time scale and formula 
(jlj) is replaced with 

o 

(14) m 



J -I 

Next, new variables are introduced so as to capture evolution up to t ~ 1/e, 
X = ex and the slow time T = et. 

Note that a solution to the eikonal equation is in general not smooth 
globally in time, unless the operator L has constant coefficients and the 
phase is linear, of the form {p(t, x) = P ■ x — ut. 

Formal three-scale asymptotics for solutions u'^ to (j2I) are given by Donnat 
l,Don94l. 



(15) u"(t,a;) ~ e'"^e"t/n(T,X,t,x 



5 ^ j\{T,X,e)=(et,ex,{l3-x-ujt)/e)- 



n>0 



The natural dependence of the profiles Un (periodic in 9) in the variable X 
(and x) is of Sobolev type, and continuous (thus bounded) in T G [0,To], 
whereas for consistency of the Ansatz, i.e. \eUn+i\ -C \Un\-, suhlinear growth 
in t = T/e is required, 

Vn>l, l||f/„||^^0. 

Donnat, Joly, Metivier and Ranch jDJMR96] give a rigorous justification 
of this monophase approximation, for semilinear hyperbolic (non-dispersive) 
systems, and profiles with no non-oscillating part (assuming that the Taylor 
expansion of the nonlinearities contains odd parts only). The Schrodinger 
profile equation (see (fT5|) ) is interpreted as a diffractive correction to geo- 
metric optics on long times, by analogy with Fresnel's diffraction. 

In order to get faster numerics for nonlinear optics, Colin, Gallice and 
Laurioux (Barrailh) |CGL05j introduce intermediate models between geo- 
metric and diffractive optics, keeping a dependence of the eigenvalue A on 
£, instead of Taylor expanding up to the order 2, so that the group velocity 
becomes a pseudo-differential operator v{k + ed^). 
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Rectification. For general systems, even if tlie initial data is purely os- 
cillatory, a non-oscillating mode (or mean field) will be created. This is 
optical rectification. In |JMR98j , Joly, Metivier and Ranch allow such a phe- 
nomenon, still for non- dispersive systems. The case of semilinear dispersive 
systems is treated by Lannes |Lan98j . We now illustrate the difficulty caused 
by interaction of oscillating and non-oscillating modes, and simply look at 
the obtention of the first profile equations. 

In the dispersive case, ttq 7^ Id (the origin belongs to the characteristic 
variety C of L) and the set of nonzero eikonal frequencies, 

S = {aeZ* \ a{-uj,k) G C}, 

is finite, in general. We assume that the only singular point of C in {a{—uj, k) \ 
a E S} is the origin. Splitting profiles Un = X^aez (^5 ^5 '^j ^) ^^""^ i^^o 
average and oscillations U* = Un — U^, one gets from the first powers of 
e in the formal expansion of Lu^ — F{u'^), 

(16a) n = Ul vro f/° = Ul 

(16b) 7roLi(9t,,)7ro t/° = 0, 

(16c) 7roLi(9T,x)vro U^^ - 7roLi{dt,.)L^'Li{dt,,)7ro 

+ 7^oFo(f/o^ f/o°) = -{dt + ■ d^) f/f , 

and Va G Z\{0}, 

(17a) (dt + v^-d,) U^ = 0, 

(17b) {dr + Va ■ dx) K - nMdL.)Lita(3)-'L,{dt,,)n^ 

+ -nMUl U^) = -{dt + Va ■ d,) f/f . 

As for the transport operator at the group velocity, an "algebraic lemma" 
|Tex04j shows that the second order operator HaLi{dt,x)L{iaP)~^Li{dt,x)T^a is 
in fact |9|A(a/3) • {d^, dx). On the other hand, the operator TiQLi{dt^x)T^Qi cor- 
responding to low frequencies, and called "long-wave operator" , is symmetric 
hyperbolic, but is not a transport operator, since the origin is a singular point 
of the characteristic variety C of L. We simply know |Lan98j that its char- 
acteristic variety Clw is the tangent cone to C at the origin. This cone may 
contain hyperplanes {(r, ^) G M^^"^ | t + ■ i = 0}; this is precisely the 
rectification criterium. 

The difficulty now lies in the compatibility of the profile equations. Lannes 
|Lan98j defines nice average operators, which are the analytic analogue of the 
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algebraic projector 11. They identify and separate the various propagation 
modes at scale (T, X) and give necessary and sufficient conditions for the 
profile equations to be solvable, with a t-sublinear corrector Ui. Split the set 
of eikonal frequencies S into the resonant set 

SR = {aeS\{T + v^-^ = 0} C Clw}, 

-denoting f/g'" the corresponding modes of Uq, and Uq'" , a' G the others- 
and the non-resonant set, Snr = £\£r. Now, the abovementioned conditions 
consist in replacing (jl6c|) . ()17b|) . for a resonant mode a G £r, with 

(18a) 

7roLi(c'T,x)vro ^/q'" - 7roLi(9i,,.)Lo ^Li(9i,,.)7ro f/°'" + 7roFo(f/o°, t/o'") = 0, 
(18b) 



{Ot + ■ dx) - -9|(a/3) ■ (9,., d,) + n^K^ f^o'°) = 0, 

whereas equations for non-resonant modes (a G Snr, ct' ^ are decoupled 
from the others, 



(19a) 

7roLi((9r,x)7ro U^'"' - 7roLi(9t,x-)^o ^^i(^t,x)7i"of^o'°' + '^oFo{0, Uq'"') = 0, 
(19b) 

{dr + ■ dx) - '-dl{af3) ■ [d,, d,)U^ + 7r^F^{U^, 0) = 0. 

The corrector UUi absorbs the difference between these equations and (jl6cp . 
(IlTbl), 



{dt + ■ d,) f/f = n^F^iU^, U^) - F^{U^, f/°'")], Va G S. 

This mode decoupling produces a solvable system of profile equations 
(polarization ()16a|) . linear hyperbolic evolution ()16b|) . ()17ap . slow nonlinear 
evolution (fTH|) ). (fTI?|) '). and clearly explains the rectification effect: even if Uq 
vanishes initially, non-oscillating modes may be created by nonlinear inter- 
actions in (jl8aj) . The asymptotics is then only at first order: the profile Uq 
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is constructed, together with correctors Ui, U2, whose secular growth (see 
|Lan03bp prevents from higher order asymptotics. 

The analogue study for hyperbolic systems with variable coefficients is 
done by Dumas |Dum03b] . |Dum04j . for WKB expansions with several non- 
linear phases: a rescaling in ()13p leads to a "weakly nonplanar" version of 
long-time diffraction, 

u'^{t,x) ~ e"^ e^Uniet, ex, ^{et, ex) / e , 0{et, ex) / e^) . 

n>0 

The validation of such asymptotics necessitates coherence assumptions on the 
phases 0, ip relatively to the operator L as well as to the tangent operators 

TTaLTCa- 

Self-focusing. Solutions to a (focusing) nonlinear Schrodinger equation 
may blow up in finite time, and this is usually interpreted in laser physics as 
the self- focusing of the laser beam jSS99j . Now, if the solution to the original 
Mawxell equations describing the beam is globally defined, what does this 
singularity mean? In |Dum03aj . Dumas shows that the diffractive profile 
blow-up corresponds to focusing for a perturbed weakly nonlinear geometric 
optics model, and Schrodinger approximation is valid at least up to a time 
of the order of a negative power of ln(l/£:) before blow-up. This shows that, 
even if each solution is defined globally in time, it undergoes, between 
t = and t = ts, an amplification by a positive power of \ia{l/e). 

Pulses and continuous spectra. The propagation of pulses over diffrac- 
tive times is described by Alterman and Ranch in the linear |AR00j and 
nonlinear case jAR03j . The profiles from (jlSp then have compact support in 
9. An important difference with diffractive optics for oscillatory wavetrains 
resides in the profile equations 

drd^U - ^d^Xmidl + dl)\J = F{U). 

The operator is not defined on the space of profiles considered, which can 
be seen, by Fourier transform, as a small divisors problem. Alterman and 
Ranch solve this difficulty using "infrared cut-offs" : they define approximate 
profiles by truncating the low frequencies in the equation, and then show 
the convergence of the approximate solution based on these approximate 
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profiles. See also |SU03j : even if Schneider and Uecker do not compare 
their results with the ones from the WKB method, they address the problem 
of existence and stability of diffractive pulse solutions to nonlinear optics 
Maxwell's equations, using a center manifold reduction (so that they get 
exponential asymptotic stability, but only around some particular family of 
solutions) . 

In |BL()2j . Barrailh and Lannes extend this approach to profiles with 
continuous spectra, which model Raman scattering and lasers with large 
spectrum. 

6.3 Transparency and larger amplitudes 

Transparency. The analysis above enlights the nonlinear interactions lead- 
ing to rectification. Unfortunately, computations on physical models, such as 
Maxwell-Bloch equations , or ferromagnetic Maxwell equations, reveal that 
the nonlinearities in (jl8a|) simply vanish! This phenomenon, called weak 
transparency, is expressed as. 



In order to reach nonlinear regimes, one may then increase the observation 
time, or the amplitude. The latter is studied (for geometric optics time C(l)) 
by Joly, Metivier, Ranch in [jMROObJ for semilinear systems of Maxwell- 
Bloch type (including the physical ferromagnetism system in space dimension 
one). 



with L and M symmetric hyperbolic operators as in Q, and /, g, h bilin- 
ear nonlinearities. When weak transparency occurs, they look for solutions 
with larger amplitude than the usual ones {0{1) instead of e). They show 
that weak transparency is necessary to construct formal WKB expansions, 
and give a necessary and sufficient criterium (strong transparency) for the 
stability of these WKB solutions. Strong transparency and the particular 
structure of the system in fact allow a nonlinear change of unknowns which 
brings back to the weakly nonlinear setting. 



V?7eC^,VaeZ, IT, 



a 



5^F,(7r^?7,7r,_^f/) =0. 



(20) 
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The same transparency property allows Jeanne |Jea02j to construct geo- 
metric optics asymptotics of large solutions to (semilinear) Yang-Mills equa- 
tions from general relativity. 

Using Joly-Metivier- Ranch's method, Colin |Col02j rigorously derives Davey- 
Stewartson (DS) systems (coupling of (NLS) and a hyperbolic or elliptic 
equation, modeling for example shallow-water, starting from Euler equa- 
tions with free surface) from Maxwell-Bloch type systems, over diffractive 
times. In space dimension one, Schneider |Sch98aj . |Sch98bj obtains similar 
results via normal form techniques. Colin also shows that the obstruction for 
such a derivation from general hyperbolic systems comes from rectification 
effects. Next, Colin and Lannes |CL04j perform the extension to general sys- 
tems, and apply their results to the Maxwell-Landau-Lifshitz ferromagnetism 
model (see Landau and Lifshitz jLL69j), getting mean- field generation. This 
corresponds to the physics papers of Leblond on pulse propagation, deriv- 
ing (NLS) [LebOlj . combining with an expansion of optical susceptibilities 
Leb02j . or deriving (DS) in ferromagnetic media |Leb96j . jLeb99j . 

For ill-prepared data allowing rectification, they build a (DS) approxima- 
tion (valid only over times of the order of ln{l/e)). To this end, they need a 
long- wave correction to the original Ansatz. Furthermore, they assume (in or- 
der for the limit (DS) system to be well-posed, i. e. to be a Schrodinger-elliptic 
coupling) that long-wave-short-wave resonance (CROLOC, in French) does 
not occur, which means that the characteristic varieties of the long-wave op- 
erator and of the original operator are not tangent away from the origin. In 
order to reach 0{l/e) times with rectification, they consider |CL01j solutions 
with size y/e (intermediate between 1 and e), and use four-scale profiles to ob- 
tain their CROLOC. A numerical study of long-wave-short-wave resonance 
is due to Besse and Lannes |BL01j . 

Some works on waves in plasmas rely on similar techniques. The basic 
model is then the quasilinear Euler-Maxwell system (see Sulem and Sulem 
|.SS99] ). Cohn, Ebrard, Gallice and Texier |CEGT04] study a simpler model, 
a Klein-Gordon-wave coupling, for which a change of unknowns leads to a 
semilinear system with weak transparency property. The diffractive time 
approximation is the Zakharov system, for electromagnetic field u, and ion 
population n, 

idfU + Au = nu, 
din — An = A|Mp. 

For the full Euler-Maxwell system, Texier |Texn5j rigorously derives, in the 
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geometric optics regime, a weak form of the above system (where e stands 
in front of the terms Am and nu, thanks to "generahzed WKB asymptotics" , 
i.e. without ehminating the residual from profile equations. 



Conservation laws with a linearly degenerate field. Cheverry, Gues 
and Metivier have adapted the notion of transparency above to quasilinear 
systems of conservation laws. They classify nonlinear regimes on times 0{1) 
with asymptotics 

x) ~ uo(t, e^'^Ukit, X, LP Is). 

k>l 

The value / = 1 corresponds to weakly nonlinear regime; / = 2, to strong 
oscillations] Z > 3, to turbulent oscillations; / = cxd, to large amplitude 
[i.e. 0{1)) oscillations (more singular solutions are the stratified solutions of 
Ranch and Reed |RR88j . studied in the quasilinear case by Corli and Gues 
|GG01p . A difficulty in WKB analysis, when dealing with amplitudes larger 
than the weakly nonlinear ones, comes from the fact that the hierarchy of 
profile equations changes. In particular, a coupling appears between phase 
and amplitude in the eikonal equation. 

In |CGM03j . they show that linear degeneracy of one of the fields is a 
weak transparency condition that ensures existence of formal strong solu- 
tions. In [ GGM04j . they study existence and stability of large amphtude 
(polarized) waves, under a stronger transparency assumption (existence of a 
good symmetrizer and constant multiplicity of the linearly degenerate eigen- 
value) that relates the system to Euler non-isentropic gas dynamics equations 
-see Metivier and Schochet |MS01 I| about the stability of large oscillations 



in time for the entropy. Earlier results on this topic were only formal ones 
(Serre |Ser95j ). or in space dimension one: see Peng |Pen92j . Heibig |,Hei93j . 

Corli and Gues [CX^Olj . and Museux |Mus()4j . 

Cheverry continues this work, in connection with turbulence. In |Che04j . 
he proves that (vanishing) viscosity can compensate the lack of transparency. 
In jChear j . he converts the stability problem into a "cascade of phases" phe- 
nomenon, where the phase has an asymptotic expansion, whose coefficients 
are coupled with amplitudes. 

Supercritical WKB solutions to (NLS). Coupling between phase and 
amplitude also occurs for large amplitude oscillating solutions to nonlinear 
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Schrodinger equations. Existence of such solutions before caustics is obtained 
by Grenier |Gre98j for a single nonlinear phase with infinitely accurate 
WKB expansions (extending results of Gerard jGer93j ). thanks to a "per- 
turbed phase" technique, ~ Tlk^^'^k- 

7 Long waves 

We finally mention the situation a priori opposite to highly oscillatory prob- 
lems, where wavelength is large. This is the typical framework of water waves, 
modelled by Euler equations with free surface, (EPS); see |SS99j . However, 
here again come into play the coupling with a mean field, and the "long wave 
operator" , as for diffractive optics rectification effects from previous section. 
This explains, for example, the formation of pairs of waves travelling in op- 
posite directions. 

Due to symmetries, the second-order differential terms vanish, in long- 
time asymptotic models, and the limit equations take the form of Korteweg- 
de Vries (KdV) equations. For one- dimensional quadratic quasilinear disper- 
sive systems of type ((21) , the appropriate Ansatz is 

n>0 

and profiles satisfy transport equation at scale so that, 

UQ{eH,t,x) = Vo{T,y)iT=eH,y=x±t- 

Furthermore, they are solutions to (KdV) in variables {T,y), 

drVo + ^d^Vo + ^dyiV,') = 0. 

Ben Youssef and Colin |BY('00j . as well as Schneider |Sch98cj . derive this 
asymptotics for simplified general systems, getting decoupled equations for 
each mode. The case of (EES) is treated by Schneider and Wayne with 
SW02j or without _SW00j surface tension. The case of transverse pertur- 
bations is more singular, since it leads (for profiles U{eH,ex2,t,xi)) to a 
Kadomtsev-Petviashvili-type equation (KP), 

drv + ^d^v + \d-^dlv + \dy{v') = 0. 
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This is obtained by Gallay and Schneider for unidirectional waves, and then 
by Ben Youssef and Lannes |BYL02j for pairs of waves. Note that (KP) 
equations are singular with respect to low frequency (in y) solutions. The 
"infrared cut-off" technique from section 16.21 is thus useful here -and this 
singularity poses the problem of consistency of the approximate solution, 
since residuals Lm^ — F{u'^) may not be small (see Lannes |Lan03aj ). 

The equations are asymptotically decoupled, but coupling appears at least 
when propagation takes place in a bounded domain (or for periodic solu- 
tions). Taking coupling into account also improves convergence rates. Ben 
Youssef and Lannes highlight such coupling effects between (KP) equations 
simply using "generalized WKB expansions" as described in paragraph 16.21 
A well-known coupled approximation of (EFS) is the Boussinesq system. In 
fact. Bona, Colin and Lannes |BCL03j obtain, in 2 or 3 space dimension, a 
result ensuring the simultaneous validity (or non-validity) of a whole three- 
parameters family of such systems, including Boussinesq. For numerics on 
this topic, see Labbe and Paumond LP04J. 

The same regime is of interest in micromagnetism, i.e. for Maxwell- 
Landau-Lifshitz system (where the medium responds to the electromagnetic 
field via magnetization). Colin, Galusinski and Kaper |CGK02j study the 
propagation of pairs of travelhng waves in space dimension one, and de- 
rive a semilinear heat equation, whereas for two-dimensional perturbations, 
Sanchez |San05j obtains Khokhlov-Zabolotskaya equations. 
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